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A brief non-technical review of the recent study Q of classical integrable structures in 
quantum integrable systems is given. It is explained how to identify the standard objects of 
quantum integrable systems (transfer matrices, Baxter's Q-operators, etc) with elements of 
classical non-linear integrable difference equations (r-functions, Baker-Akhiezer functions, 
etc). The nested Bethe ansatz equations for ^fc_i-type models emerge as discrete time 
equations of motion for zeros of classical r-functions and Baker-Akhiezer functions. The 
connection with discrete time Ruijsenaars-Schneider system of particles is discussed. 



At present it is argued that classical and quantum integrable models have much deeper interrela- 
tions than any kind of a naive "classical limit". In [|^, a particular aspect of this phenomenon has 
been analysed in detail: Bethe equations, which are usually considered as a tool inherent to quantum 
integrability, arise naturally as a result of solving entirely classical non-linear discrete time integrable 
equations. This suggests an intriguing connection between integrable quantum field theories and 
classical soliton equations in discrete time. 

Here we outline main ideas of the paper Q omitting all technical details. To save the space, the 
material is organized as a "quantum - classical dictionary" supplied with brief comments. 

1. Eigenvalues of quantum transfer matrices = (classical) r-functions. Due to the Yang- 
Baxter equation the transfer matrices commute for all values of the spectral parameters in the auxiliary 
space (AS)|]: [r4(M), T4'(m')] = 0- This property allows one to diagonalize them simultaneously. 
From now on we use this diagonal representation. The identification with r-function is justified in the 
J> .next item. 

. 2. Fusion rules = Hirota's difference equation. The fusion procedure in the AS gives rise to 
5^ I a family of commuting transfer matrices Tji^{u) with the same quantum space. They obey a number 
of fusion relations |Q] which can be recast into the model- independent bilinear form |^]. Let T"(n) be 
the transfer matrix for the rectangular Young diagram of length s and height o, then it holds 

T^{u + l)T^{u - 1) - Tt^,{u)T^_,{u) = T^-'\u)Tr\u) . (1) 

Remarkably, this equation coincides with Hirota's bilinear difference equation (HBDE) |^ which is 
known to unify the majority of soliton equations, both discrete and continuous. 

Fusion of more complicated representations in the AS is described by higher representatives of 
the hierarchy of HBDE-like equations 0. For example, consider Young diagrams consisting of two 
rectangular blocks (i.e. with ai lines of length si + S2 and the rest 02 lines of length si) and let 
Tg^^f^{u) be the corresponding transfer matrix. Then it holds 

+ K'^i:,' (u - i)t:^^\]':^- [u + 1) 

= c+m"/ - + 1) + ci? - + 1) • (2) 



^Talk presented at the II Sakharov International Conference 

■^Here and below we use the notions and terminology of the quantum inverse scattering method, see e.g. 
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In what follows we consider eq. ([T|) only. 



3. Specifying a particular quantum model = imposing particular boundary and ana- 
lytic conditions in HBDE. 

For models associated to A^-i-type quantum algebras these conditions are: 

Tg{u) = as a < or a > fc; 

T^(u) = for any - A; < s < 0, < a < A; . (3) 

It follows from (Q), ^ that Tg{u) and Tg{u) factorize into a product of functions of n + s and u — s. 
We adopt the normalization in which Tg{u) = (/)(n + s), Tg{u) = (f){u — s — k), where the function <i){u) 
carries information about the particular quantum model. 

Another important condition (which follows, eventually, from the Yang-Baxter equation) is that 
Tg(u) for models with elliptic i?-matrices has to be an elliptic polynomial in the spectral parameter u. 
(By elliptic polynomial we mean essentially a finite product of Weierstrass o"- functions.) For models 
with rational i?-matrix it degenerates to a usual polynomial in u. 

4. "Step" of nested Bethe ansatz = Backlund transformation. Quantum integrable models 
with internal degrees of freedom can be solved by the nested (hierarchical) Bethe ansatz method. The 
method consists essentially in integration over a part of degrees of freedom by an ansatz of Bethe type, 
the effective hamiltonian being again integrable. Repeating this step several times, one reduces the 
model to an integrable model without internal degrees of freedom which is solved by the usual Bethe 
ansatz. 

The classical face of this scheme is Backlund transformation, i.e. passing from solutions to the 
non- linear equation to (properly normalized) solutions of the auxiliary linear problems (ALP), which 
satisfy the same non-linear equation. The ALP for HBDE have the form Q 

T^^I{u)F''{s, u) - T^+\u - l)F'^(s + 1, n + 1) = T^{u)F''+\s + 1, n) , 

riVi(n - l)F-{s,u) - Tt{u)F\s + l,u-l) = T^+\u - l)F--\s + l,u), (4) 

Eq. (HI) is the compatibility condition for (^). Due to the symmetry between T and F [|8|, the latter 
satisfies the same non-linear equation (|l|), so the transition T ^ F is a Backlund transformation. The 
important point is that the boundary condition for F"'{s,u) is the same as in (^), the only change 
being a reduction of the Dynkin graph: /c — > — 1. In other words, the number of non-zero functions 
gets reduced by 1. Using this property, one can successively reduce the j4fc_i-problem up to Ai. 

5. Eigenvalues of Baxter's Q-operators = properly normalized solutions to the ALP 
(Baker- Akhiezer functions). To elaborate the chain of Backlund transformations, let t = 0, 1, . . . ,k 
mark steps of the flow Ak_i — > Ai and let F^_^^{s,u) be a solution to the ALP at (/c — t)-th step (in 
this notation F^{s,u) = Tg{u)) such that F^{s,u) = as a < or a > t. Due to this condition F^ 
and Fl are "chiral" functions, i.e. 

F^{s,u)=Qt{u + s), F^{s,u) = Qt{u-s-t), (5) 



where Qt (u) are some functions playing a distinguished role in what follows since they can be identified 
with generalized Baxter's Q-operators in the diagonal representation (see below). Here is an example 
of this array of r-functions for the A2-case (A; = 3): 

1 

Qi{u + s) Qi(u-s-l) 

Q2iu + s) F^{s,u) Q2{u-s-2) 

Q-iiu + s) F^{s,u) F3is,u) Q3{u-s-3) 
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(in this case Qz{u) = (j){u)). 



6. Nested Bethe ansatz equations = Calogero-type models in discrete time. It follows 
from that t^'°'{u) = F^_^{u + a,u) satisfies the bilinear equation 



r 



*+1''^(u)t*'"+^(u) - T*'"(n)r*+^''^+i(n) = T*+^'''(n + 1)t*''^+^(u - 1) (7) 



which is HBDE in "light cone" variables. At the same time t^'^{u) = Qk-tC^u). Let u*- be zeros of 
Qt{u): Qt{u'j) = 0. For models with elliptic i?-matrices Qt{u) should be elliptic polynomials in u. This 
condition leads to a number of constraints for u*- which can be derived using the technique developed 
in 1^] for elliptic solutions to the KP equation. These constraints are nothing else than the nested 
Bethe ansatz equations: 

Q,-i(^ + 2)Q,(^-2)Q,+i(4) ^ _ 

Qt^i{u])Qt{u'j + 2)Qt+i{u]-2) ^' 

(with the boundary condition Qo{u) = 1, Qk{u) = 4>{u))- These equations can be understood as 
"equations of motions" for zeros of Qt{u) in discrete time t (level of the Bethe ansatz which runs over 
the Dynkin graph). The analogy with elliptic solutions of the KP equation suggests to call them the 
discrete time analogue of the Ruijsenaars-Schneider (RS) system of particles (see also |]lO|). Taking 
the continuum limit in t (provided the number of zeros Mf = M of Qt{u) in a fundamental domain 
does not depend on t), one can verify that this system of equations does yield equations of motion for 



the RS system |11| with M particles. 

However, integrable systems of particles in discrete time have a richer structure than their con- 
tinuous counterparts. In particular, the total number of particles may depend on the discrete time. 
Such a phenomenon is possible in continuous time models only for singular solutions, when particles 
can move to infinity or merge to another within a finite period of time. Remarkably, this appears to 
be the case for the solutions to eq. ^ corresponding to eigenstates of quantum models. It is known 
that the number of excitations at t-th level of the nested Bethe ansatz solution does depend on t. In 
other words, the number of "particles" in the associated discrete time RS system is not conserved. At 
the same time the numbers Mj may not be arbitrary. It can be shown that for models with elliptic 
i?-matrices in case of general position Mt = {N/k)t, where is the number of sites of the lattice 
(degree of the elliptic polynomial <l){u)). In trigonometric and rational cases the conditions on Mt 
become less restrictive but still they may not be equal to each other. 

At last we should indicate how to identify our Qj's with Qt^s from the usual nested Bethe ansatz 
solution. This is achieved by the following factorization formula which allows one to express Tg{u) in 
terms of Qts: 



E(-i 



a=0 
29„ 



,a-k Tliu + a-l) 2ad^ 

(l){u - 2) 



Qk{u)Qk-i{u-2)\ ( ^2du Q2{u)Qi{u-2) \{^Q^ Qi{u) 



Qk{u-2)Qk-i{u)J"\ Q2{u-2)Qi{u)J\ Qi{u 

where the shift operator e^" acts as usual: e^'^f{u) = f{u + 1). For the proof see ref. This formula 
coincides with the one known in the literature (see e.g. ||12[). The left hand side is known as the 
generating function for Tf (n). 

There is no doubt that this dictionary can be continued to include not only spectral properties of 
quantum system on finite lattices (which just correspond to elliptic solutions of HBDE) but dynamical 
properties and correlation functions as well. Perhaps this amounts to considering solutions to the same 
Hirota equation of a more complicated type. Another intriguing problem is to extend this dictionary 
to off-diagonal elements of quantum monodromy matrices and quantum L-operators and identify them 
with some objects in hierarchies of classical integrable equations. 
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